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Abstract 

The Hamiltonian eigenvalue problem for bound states of a quantum field theory is formulated 
in terms of Dirac's light-front coordinates and then approximated by the exponential-operator 
technique of the standard coupled-cluster method. This approximation eliminates any need for the 
usual approximation of Fock-space truncation. Instead, the exponential operator is truncated and 
the terms retained are determined by a set of nonlinear integral equations. These equations are 
solved simultaneously with an effective eigenvalue problem in the valence sector, where the number 
of constituents is small. Matrix elements can be calculated, with extensions of techniques from 
standard coupled-cluster theory. 



a Presented at the Sixth International Conference on Quarks and Nuclear Physics, April 16-20, 2012, Ecole 
Poly technique, Palaiseau, Paris. To appear in the proceedings. 
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I. INTRODUCTION 



The light-front coupled-cluster (LFCC) method [l| is intended as a method for the calcu- 
lation of hadron structure in terms of Fock-state wave functions without the usual truncation 
of Fock space. The Fock-state representation of a hadron, such as a proton, takes the form 

\p) = i) U ud\uud) + i) U udg\uudg) + i) U ud gg \uudgg) + ip uudq q\uudqq) H . 

This is to be an eigenstate of the QCD Hamiltonian 

(K.E. + Vq CD ) \p) = E p \p) 

with E p = y/rrip + p 2 and Vqcd the interaction terms of gluon emission and absorption by 
quarks, quark-antiquark production and annihilation, and the three and four-gluon vertices. 
This is equivalent to a set of coupled integral equations for the wave functions. 

Usually, this infinite set of equations is made finite by truncation of Fock space; however, 
this introduces uncanceled divergences [2j. For example, the Ward identity of gauge theories 
is destroyed by truncation, because the limitation on the number of particles allowed in 
intermediate states removes some of the contributions to the identity. There is a direct 
analog in Feynman perturbation theory, where separation of a covariant diagram into time- 
ordered diagrams, and removal of those that include intermediate states with more particles 
than some finite limit, destroys covariance, disrupts regularization, and induces spectator 
dependence for subdiagrams. In the nonperturbative case, this happens not just to some 
finite order in the coupling but to all orders. 

The method is formulated in terms of light-cone coordinates 0, 0] • The time coordinate 
is x + =t + z, and the spatial coordinates are x = (x~ , x±), with x~ =t — z and x± = (x, y). 
The light-cone energy is p~ = E — p z , and the momentum is p = {p + ,p±), with p + = E + p z 

and p± = (p x ,Py). This leaves the mass-shell condition p 2 = m 2 as p~ = m ^ p± and the 
mass eigenvalue problem as jEf 



M 2 -I- P 2 

v-\p)= p + H p), v\p) = p\p). 

The advantages of this coordinate choice include the absence of spurious vacuum contri- 
butions to eigenstates and a boost-invariant separation of internal and external momenta. 
Vacuum contributions are suppressed because p + is positive for all particles; the vacuum 
cannot produce particles and conserve p + . The separation of internal momenta is obtained 
by defining longitudinal momentum fractions X{ = pj / P + and relative transverse momenta 

h± = Pi± - XiP-L- 

The original coupled-cluster (CC) method was developed for the many-body Schrodinger 
equation in nuclear physics @ and later applied to the many-electron problem in molecules 0] • 
The method has become an important tool in nuclear physics and physical chemistry for 
the iV-body problem in nonrelativistic quantum systems [8j. The basic idea is to form 
an eigenstate as e T |0), where \4>) is a product of single-particle states and the terms in 
T annihilate states in \(f>) and create excited states, to build in correlations; however, the 
number of particles does not change. The approximation made is to truncate T at some 
number of excitations. The LFCC method uses the mathematics of these constructions but 
applies it to a situation where \<p) contains a small number of particles, T adds additional 
particles, and the states are eigenstates of momentum with Dirac-delta normalization. 
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II. LIGHT-FRONT COUPLED-CLUSTER METHOD 



To solve the fundamental eigenvalue problem V~\ip) = — p+ x \tp) by the LFCC method 
we write the eigenstate as \ip(P)) = \/~Ze T \4>{E)) an d seek solutions for the valence 
state \<f)(P_)) and the operator T. This operator contains terms that only increase par- 
ticle number, while conserving J z , light-front momentum P, charge, and the other quan- 
tum numbers of the eigenstate. The constant Z controls normalization, which is chosen 
to be (ip(P!)\ip(P)) = 6 (Pi — P), with the valence state normalized in the same way: 
(<fi(P')\<f>(P)) = S(P? — P). Because p + is positive, T must include annihilation, and powers 
of T include contractions. This converts the original eigenvalue problem into an eigenvalue 
problem for the valence state 

W- _|_ P 2 

p v v-\m) = p + ± \m)), 



with V— = e T V e T the effective Hamiltonian and P v the projection onto the valence Fock 
sector, and into an auxiliary equation for T 

(l-P„)Pi </>(£)> = 0. 

Calculation of expectation values, and more generally matrix elements, requires some 
care, to avoid any necessity of computing the infinite sum implied by the inner product 
(0|e Tt e T |0). This can be done, with use of constructions from the CC method (8j. The 
expectation value 

<d)\e T 'de T \d)) 



(O) 



|e Tt e T l 



for an operator O is rewritten as (O) = (^\0\(j)) in terms of an effective operator O = e T Oe T 
and a left eigenvector of V~ : 



Tt T 

e e 



|e Tt e T l 



The effective operator can be computed from the Baker-Hausdorff expansion = + 
[O , T] + |[[0, T], T] H — • . The bra (ip\ is seen to be a left eigenstate by the following steps: 



e Tt P"e T — t e Tt e T M 2 + P]_ , 7l 

P_ ( "' ' '\e T 'e T \6) = (<b\e T 'e T \6) = ~ 



Also, it has the normalization 

${Z)\m) = ^)\ ,/ T f T ux \^p)) = ^e-p). 

As formulated to this point, the new valence eigenvalue problem and the auxiliary con- 
dition for T provide an exact representation of the original eigenvalue problem. The LFCC 
method then invokes a truncation, not of Fock space but of the terms in T and in the pro- 
jection 1 — P v . This leads to a finite set of auxiliary equations for a finite set of functions 
that define the truncated T operator and to a finite number of terms in the valence pro- 
jection of the effective Hamiltonian V~. The latter is conveniently expressed in terms of a 
Baker-Hausdorff expansion where the number of contributing terms is finite. Similarly, the 
contributions to matrix elements of operators are also finite in number. 
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III. APPLICATION TO A SOLUBLE MODEL 



To illustrate the method, we apply it [1] to an exactly soluble model |9j, a light-front 
analog of the Greenberg-Schweber model with a static fermionic source that emits and 
absorbs bosons without changing its spin. The model lacks full covariance because of the 
static source; all states are limited to having a fixed total transverse momentum P±, which 
we take to be zero. In this context, not all features of the method are apparent, but the 
model is sufficient to show how the method is used. 

The light-front Hamiltonian is [l[ 

V = [ dp M2 + ^ P+ J2bi(p)b s (p) + / dg^y^^aj^ig) (3.1) 

J s J I Q 

+ ^T&E" G^f) 7 [ a Kin(gMp + g) + bt(p + q)b s (p)a l (g)_ , 

where aj creates a "physical" boson of mass /io, a\ creates a Pauli-Villars (PV) boson of mass 
fii, and b\ creates the fermion with mass M and spin s. The valence state is the bare-fermion 
state \4> a {P_)) = &J.(P)|0). The T operator is truncated to a single boson emission: 

T = 5Z/ dqdpt ls (q,p)al(q)bl(p)b s (p + q), 

Is 

and the projection 1 — P v is truncated to the one-fermion/one-boson sector. 
This form generates the exact solution, with 

The fermion self-energy contribution M' is the same in all Fock sectors and the effective 
Hamiltonian |l| contains all three of the contributions analogous to those for the Ward 
identity in QED. The fact that the self-energy loop is the same in the valence sector and 
the one-fermion/one-boson sector plays a critical role in yielding the exact solution. 

To compute an observable, we consider the Dirac form factor for the dressed fermion. It 
can be obtained from a matrix element of the current J + = ipj+tp coupled to a photon of 
momentum q. The matrix element is generally [Io[ 



with Fi and F 2 the Dirac and Pauli form factors. In the model, the fermion cannot flip its 
spin; therefore, F 2 is zero, and we compute only F\. Also, there are no contributions from 
fermion-antifermion pairs, so that the current is simply 

In the LFCC method, the form factor is approximated by the matrix element 



F, (q 2 ) = Svr 3 (P + q) \ J+ (0) 1 ± (P) ) , 
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with J + (0) = J + (0) + [J + (0), T] + ■ ■ ■ . The truncated left-hand eigenvector is 

(r(p)\ = (r(p)\ + E/ - « + )Cfe£)<ok(9)6.(z- ?), 

Zs 

where is the left-hand one-fermion/one-boson wave function. If this wave function is 
assumed to take the form 

substitution into the left-hand eigenvalue problem yields a one-dimensional integral equation 
for l(y) 

Kv) = i + #2 4^ f dy'(i - yO - y) 2 W - *)) - %')]• 

1D7T /io/^l JO 

The solution of this equation can then be used to compute matrix elements. 

Only the first two terms of the Baker-Hausdorff expansion of J + {0) contribute to the 
matrix element. The first term contributes 1/87T 3 and second contributes 

(^ ± (Z + g)|[J + (O),r]|0 ± (P)) = -^^(-iy I dq>6(P + + q + - q' + ) 

x^±W,P + q)[0(P + ~ q' + )ti±W,P- i) - t l± tf,P + q - g')l- (3-2) 
The form factor is then 

{ 2 I" rl/(l+a) 
16tt 2 ^ ' ~' fjfifil J 

- [ dyi(y)y(l-y)*> 
Jo 



dyl(y)y(l-yy[l-(l + a) y y 



(3.3) 



with a = q + / P + . In the limit of q 2 — > 0, we have a = and -Fi(O) = 1, which is exactly the 
unit charge in the current J + = 



IV. SUMMARY 

The advantages of the LFCC method are the absence of Fock-space truncations that cause 
uncanceled divergences, elimination of Fock-sector and spectator dependence of self-energy 
contributions, and provision for systematic improvement, through the addition of terms to 
the truncated operator T. The terms in T can be organized according to the nature and 
number of particles annihilated and created. Applications to theories beyond the simple 
model considered here are in progress, with some preliminary work on the dressed-electron 
state in QED already completed Additional work in QED will include consideration 
of the dressed-photon state, extension of the dressed-electron state to include e + -e~ pairs, 
muonium, and positronium. Application to QCD will begin with consideration of mesons in 
holographic QCD ' 
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